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Abstract. We estimate the extragalactic muon neutrino background which is produced by hadrons injected by 
very young pulsars at an early phase after supernova explosion. It is assumed that hadrons are accelerated in 
the pulsar wind zone which is filled with thermal photons captured below the expanding supernova envelope. In 
collisions with those thermal photons hadrons produce pions which decay into muon neutrinos. At a later time, 
muon neutrinos are also produced by the hadrons in collisions with matter of the expanding envelope. We show 
that extragalactic neutrino background predicted by such a model should be detectable by the planned 1 km 2 
neutrino detector if a significant part of pulsars is born with periods shorter than ~ 10 ms. Since such population 
of pulsars is postulated by the recent models of production of extremely high energy cosmic rays, detection of 
neutrinos with predicted fluxes can be used as their observational test. 
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1. Introduction 

The construction of high energy neutrino telescopes 
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^ \ with the effective area > 10 5 m 2 (AMANDA^IceCube, 
K> ' ANTARES, NESTOR) will probably allow observation of 
rS \ first neutrinos of astrophysical origin, giving for the first 
■ time clear information about high energy hadronic pro- 
' cesses in the universe. Therefore estimations of neutrino 
fluxes from discrete sources and their contribution to the 
extragalactic neutrino background (which may be proba- 
bly easier to detect) is so important. 

Up to now the flux of extragalactic neutrino back- 
ground has been estimated with rather high confidence 
in the case of interaction of diffuse cosmic rays with the 
matter and radiation. Estimations of neutrino fluxes from 
discrete sources, where as it is well known, high energy 
processes are important i.e. active galactic nuclei and 7- 
ray bursts, are less certain. Also more exotic processes, 
like decays of topological defects, have been discussed as 
possible sources of extremely high energy neutrinos. For 
recent reviews of these estimations we refer to the articles 
by Protheroe (1998) and Learned & Mannheim (2000). 

Supernova explosions, during which young pulsars are 
formed, seem to be another likely sites of high energy 
processes in which observable neutrino fluxes are ex- 
pected (e.g. Berczinsky & Prilutsky 1978; Bednarek & 
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Protheroe 1997; Protheroe, Bednarek & Luo 1998). In our 
recent paper (Beall & Bednarek 2001, BB) we calculate 
the muon neutrino (and antineutrino) flux expected from 
a very young pulsar surrounded by the supernova enve- 
lope in our Galaxy. It has been shown that muon neutri- 
nos are likely to be observed in the large scale (1 km 2 ) 
neutrino detector from the supernova within our Galaxy 
if a pulsar with a period shorter than ~ 10 ms is formed. 
However such explosions happen relatively rarely and we 
would have to be lucky to detect such a supernova within 
reasonable period of time. Therefore in order to test the 
predictions of the mentioned model by the planned neu- 
trino detectors, we calculate in this paper the extragalactic 
neutrino flux produced by very young pulsars formed in 
the universe. 

2. Neutrinos from nuclei accelerated by pulsars 

Pulsars are produced during the type II and Ib/c super- 
nova explosions of massive stars which detailed models 
have been investigated in many papers (see e.g. Woosley 
et al. 1993). When the iron core collapses, a very hot 
protoneutron star is formed. It cools to a neutron star 
within about £ns ~ 5 — 10 s from the collapse (Burrows 
& Lattimer 1986, Wheeler et al. 2000). The rest of the 
presupernova mass is expelled with the velocity of the in- 
ner radius of the order of v\ — 3 x 10 8 cm s _1 produc- 
ing the supernova envelope. Its initial very high optical 
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depth drops fast with time. The region between the pul- 
sar and the envelope is filled with thermal radiation of 
characteristic temperature To « 3 x 10 9 K (see Fig. 8 in 
Woosley et al. 1993). This radiation is not able to escape 
through the envelope because of its high optical depth. 
Assuming that the total energy of radiation inside the en- 
velope is conserved, we can estimate the temperature of 
this radiation as a function of time during the expansion 
of the envelope T{t) — Tq [Pi / (R t + v x (t NS + i))] 3/4 , where 
Ri = 3 X 10 8 cm is the initial radius of the envelope at 
the moment of formation of the protoneutron star. When 
the pulsar forms (about 10 s after the explosion), the ra- 
diation temperature inside the envelope already drops to 
~ 5 x 10 8 K. This radiation escapes through the envelope 
when the optical depth drops below ~ 10 3 . This happens 
about ~2x 10 6 s after the explosion (see BB). 

We assume that at the very early stage the pulsar 
loses energy only on electromagnetic waves. Therefore 
its period changes according to the formula F^ s (t) = 
1.04x IO^-B^-I-Pq ms , where Po, ms Pns are the initial and 
present periods of the pulsar, and Bn is the pulsar surface 
magnetic field in units 10 12 G. We neglect the pulsar en- 
ergy losses to gravitational radiation due to the r-mode in- 
stabilities since they may become important about 1 year 
after pulsar formation (e.g. Andersson 1998, Lindblom, 
Owen, Morsink 1998). 

It is argued that pulsars can accelerate particles 
to the highest energies observed in cosmic rays in the 
slingshot mechanism (see recent papers by e.g. Blasi, 
Epstein & Olinto 2000 (BEO), De Goubeia Dal Pino & 
Lazarian 2000, Bednarek & Protheroe 2001). Following 
work by Blasi et al. (BEO), we assume that a part of the 
magnetic energy in the pulsar's wind zone accelerates iron 
nuclei. The maximum energies, which nuclei can reach in 
the wind zone, depend on the pulsar parameters 
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where n GJ = P(r LC )/(2ecP) w 3.3 x 10 n Pi 2 P m 4 is 
Goldreich & Julian (1969) density, r LC « 4.77 x 10 6 P ms 
cm, c is the velocity of light, and (3 is the coefficient de- 
scribing the efficiency of acceleration. According to the 
slingshot mechanism the acceleration of nuclei occurs very 
fast so that they do not lose energy during this stage. After 
acceleration, the nuclei move balistically through the ra- 
diation field (below the envelope) and later through the 
matter of the envelope. These particles are not frozen into 
the expanding plasma of the envelope and therefore do 
not lose energy on adiabatic process. 

We apply the spectrum of iron nuclei accelerated in 
the pulsar wind zone close to the pulsar light cylinder 
derived in Beall & Bednarek (BB). In this simple model 
the number of nuclei accelerated to energies E scales as 
a part r\ of the Goldreich & Julian (1969) density at the 
light cylinder radius. It is assumed that the pulsar with 
specific parameters injects particles within some range of 
energies which is due to the fact that the magnetic field at 
different parts of the light cylinder radius is different. The 



energy spectrum of the iron nuclei below E < Pfc, which 
injected by the pulsar at a fixed age 't' of the pulsar, is 
(see details in BB) 
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Accelerated nuclei move in the pulsar wind almost at 
rest in the wind reference frame (BEO). Therefore we can 
neglect their synchrotron energy losses. However, as we 
mentioned above nuclei interact with strong thermal radi- 
ation field suffering complete photodesintegration before 
losing significant energies on pair and pion creation. 
For the considered parameters of the pulsar and the ac- 
celeration region, nucleons from desintegration of nuclei 
lose their energy mainly on pion production. Pions decay 
into neutrinos if their decay distance ~ 7807^ cm, 
is shorter than their characteristic energy loss mean free 
path. The Lorentz factors of pions, are comparable to 
the Lorentz factors of their parent nucleons, so they move 
similarly in the pulsar wind and their synchrotron losses 
do not dominate their inverse Compton scattering losses 
(ICS) in the thermal radiation. Pions lose energy on ICS 
process mainly in the Klein-Nishina (KN) regime, but not 
very far from the Thomson regime. By comparing the pion 
mean free path for ICS losses with their decay distance 
scale Ajr, we have found that pions can decay before los- 
ing energy on ICS if the temperature of thermal radiation 
drops below T~3x 10 6 K. The temperature of radiation 
below the envelope drops to this value about idee ~ 10 4 
s after supernova explosion. At this time nuclei arc in- 
jected with the maximum energies allowed by Eq. [l] and 
cool in collisions with thermal radiation mainly on pion 
production. However when the optical depth through the 
expanding envelope drops below ~ 10 3 then the radiation 
is no more confined to the region below the envelope and 
its temperature drops much faster. Therefore we conclude 
that nucleons are able to cool efficiently in the radiation 
and produce pions which decay into muon neutrinos, only 
from <dcc ~ 10 4 s up to t con f rj 2 x 10 6 s after explosion. At 
later times, the relativistic nuclei interact with the matter 
of the envelope which density is already low enough for 
the pions to decay into high energy neutrinos. Relativistic 
nucleons find enough column density of matter in the en- 
velope (> 10 g cm -2 ) to interact at the time shorter than 
about 1 year after explosion (see BB). 

In terms of such a model we compute the spectrum 
of muon neutrinos from the pion decay, N{E V , P, P), pro- 
duced by protons accelerated by pulsars with different pe- 
riods, P, and surface magnetic fields, B (see BB). The 
pions are produced during the period 10 4 — 2 x 10 6 s after 
explosion as a result of complete cooling of nucleons in 
collisions with thermal radiation below the envelope, and 
during the period 2 x 10 6 — 3 x 10 7 s from explosion as 
a result of collisions of iron nuclei with the matter of the 
envelope. 
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3. Extragalactic neutrino background 

In order to calculate the extragalactic neutrino back- 
ground (ENB) we assume that the pulsars born at earier 
epoch do not differ significantly from those born at 
present. Thus the injection spectrum of neutrinos by a sin- 
gle pulsar depends only on its initial parameters (period, 
surface magnetic field) but not on rcdshift. Under this as- 
sumption, we can evaluate the neutrino injection rate at 
different redshifts by multiplying the injection spectrum 
of neutrinos of an average pulsar in our Galaxy by the 
neutron star formation rate, which depends on the mas- 
sive star formation rate, which in turn depends on the 
redshift. In the previous section we have calculated the in- 
jection spectrum of neutrinos from the pulsar with known 
parameters (for details see BB). Let's determine now the 
neutrino injection rate from an 'average' pulsar by inte- 
grating over the initial period and surface magnetic field 
distributions of pulsars within the Galaxy, according to 
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N(E„,P,B)f(P)—dPdB, (3) 



where P mm -r- fmax is the range of initial periods of pulsars, 
B mul and B max correspond to the range of typical values 
of the surface magnetic fields for the observed "classical" 
radio pulsars equal to 10 115 and 10 133 , respectively. 

It is believed that pulsars are born with rather short 
periods. However their exact period distribution is not 
known. Therefore we consider two possible distributions 
(see Giller & Lipski 2001). The simplest case is described 
by the uniform distribution, f(P) = l/(-P m ax — -Pmin)' And 
probably a more realistic case, described above 1 ms by 
the gamma function, which is defined by the mean pulsar 
period (P) and the index a, f(P) = AP a e- ( - a+1)p ^ p \ 
where A = ((a + 1) / (P)Y a+1 ^ / a\ after normalization to 
one. 

The distribution of surface magnetic fields of pul- 
sars at birth can be estimated from the observed pop- 
ulation of radio pulsars and assuming that the pulsar 
magnetic fields do not decay during the radio activ- 
ity stage. It can be approximated by dN/d(logB) « 
0.065-B^ year -1 , for B\2 > 2, after normalization to the 
formation rate of pulsars in the Galaxy equal to 1 per 
70 years (Narayan 1987 ). We estimate the distribution of 



pulsars with surfac e mag netic fields below 2 x 10 12 G from 
Fig. 13 in Narayan (1987). This approximation may under- 
estimate the number of pulsars with low magnetic fields, 
because of observational selection effects, and with high 
magnetic fields > 10 13 G since it does not take into ac- 
count the so called 'magnetars'. 

Now we can calculate the differential flux of ENB using 
the formula (see Fichtel & Trombka 1981), 



cN G Q(E u {l + z))F{z) 

^H o J (l + z)(l + go z)i/2 



dz, 



(4) 



where z is the redshift, Nq — 10~ 2 Mpc~ 3 is the density of 
galaxies, Q(E V (1 + z)) is the injection spectrum of neutri- 
nos with energy E v (l + z), averaged over the distribution 



of initial pulsar periods and surface magnetic fields (see 
Eq. |^). This is normalized to the pulsar birth rate per 
galaxy at the redshift z = equal to the one observed 
in our Galaxy. f(z) is the evolution factor of the pulsar 
birth rate which we take proportional to the supernova 
rate (Kobayashi, Tsujimoto & Nomoto 20001 ). Moreover 



we assume H a — 75 km s Mpc , qo — 0, and z n 

The results of calculations of ENB are shown in Figs. 1 
and 2, for the case of initial pulsar periods described by 
the uniform and gamma distributions, respectively. We in- 
vestigate how the ENB depends on the parameters which 
define those distributions. In the case of uniform distribu- 
tion of the initial periods we show the ENB from pulsars 
born within the range of periods limited by P m i n = 1 ms 
and P max = 100 ms, P m i n = 1 ms and P max = 10 ms, and 



P, 



10 ms and P n 



100 ms. In the case of gamma 



distribution we show the ENB for (p) = 10 ms and the in- 
dex a = 2, 10 ms and a = 4, and 50 ms and a = 2. In all 
cases T) = 1. Note that the ENBs, estimated for the pul- 
sars which part has been born with periods below 10 ms, 
are above the atmospheric neutrino background (ANB) 
produced in the interactions of cosmic ray particles with 
the earth's atmosphere. Therefore they can be potentially 
observed by future detectors. 

4. Discussion and Conclusion 

We considered only neutrinos produced during the time 
when the supernova envelope is optically thick. However 
pulsars inject particles also after > 1 year when the en- 
velope becomes optically thin. Without detailed model- 
ing of propagation of these particles in the supernova 
nebula it is difficult to conclude if the neutrino fluxes 
should be limited by the optically thick or thin neutrino 
bound marked in Figs. 1 and 2 (Waxman & Bahcall 1999, 
Mannheim, Protheore & Rachen 2001). The particles in- 
jected after 1 year can be captured by the magnetic field 
of the expanding large scale supernova nebula and suf- 
fer significant adiabatic energy losses (e.g. Bednarek & 
Protheore 2001). The arguments that charged particles 
convert into neutrons, e.g. protons into neutrons, (see 
Waxman & Bahcall 1999) may not be valid if the optical 
depth for hadronic collisions becomes low and adiabatic 
losses of captured hadrons are large. 

At the end of the optically thick phase and during 
the optically thin phase, the pulsar-supernova envelopes 
should become also the sources of 7-rays and neutrons. 
Production of these neutral particles has been already 
considered in a similar supernova model (with the dif- 
ference that particles are accelerated in the inner pul- 
sar magnetosphere, Protheroe et al. 1998, Bednarek & 
Protheroe 1997). The contribution of these neutral par- 
ticles to the 7-ray background and cosmic rays in terms 
of the present model will be discussed in a future paper. 

We have discussed also the problem whether the esti- 
mated extragalactic neutrino background from the whole 
population of pulsars in the universe can be detected 
by the planned experiments. The full horizontal line in 
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Fig. 1. The extragalactic neutrino background produced 
by pulsars with uniform distribution of initial periods be- 
tween 1-10 ms (dot-dashed curve), 1-100 ms (full curve), 
10-100 ms (dot-dot-dot-dashed curve) and the acceler- 
ation parameter = 1, and 1-10 ms and (3 = 0.1 
(thick dashed curve). The upper bounds on extragalac- 
tic neutrinos for the case of optically thick (MP&R, 
Mannheim et al. 2001) and thin sources (W&B, Waxman 
& Bahcall 1999) are marked by the horizontal thin 
dot-dashed lines. The atmospheric neutrino background 
(ANB) is shown by the dashed curves. The sensitivity of 
the planned IceCube detector is marked by the horizon- 



tal thin full line (Spiering pOOOj ), and the AUGER array 
by the thin dotted curves for the muon and tau v T 



neutrinos (Bertou et al. 2001) 



Figs. 1 and 2 shows the detection limit of the IceCube 
detector (Spiering 2000). It looks that ENB from pulsars 
should be easily observed by the IceCube if significant part 
of pulsars are born with initial periods shorter than ~ 10 
ms in the case of uniform distribution (see Fig. 1), and 
if the mean pulsar period is close to 10 ms in the case of 
the gamma distribution (see Fig. 2), assuming that r\ ~ 1. 
Also the planned AUGER experiment with the detection 
limits estimated by Bertou et al. 2001 may observe ENB 
from pulsars if the significant amount of pulsars is born 
with the initial periods of the order of a few ms as pos- 
tulated by the recent models for the pulsar origin of the 
highest energy cosmic rays (BEO, De Goubeia Dal Pino & 
Lazarian 2000). However, such detection would be possible 
if the hypothesis of — > v T oscillations with full mixing 
is correct. In such case the expected flux of tau neutri- 
nos, compared to the calculated flux of muon neutrinos, is 
above the sensitivity limit of the AUGER detector for the 
case of pulsars born with the milisecond periods (see lower 
dotted curve in Fig. 1, taken from Bertou et al. 2001). 



Fig. 2. As in Fig. 1 but for pulsars with the initial periods 
described above 1 ms by the gamma function with the 
mean period (p) — 10 ms and index a = 2 (dot-dashed 
curves), 10 ms and a = 4 (full curves), and 50 ms and 
a = 2 (dot-dot-dot-dashed curves) ((3=1 in the above 
cases), and for (p) = 10 ms, a = 2, and (3 = 0.1 (thick 
dashed curve). 

References 

Andersson, N. 1998, ApJ, 502, 708 
Beall, J.H., Bednarek, W. 2001, ApJ, submitted (BB) 
Bednarek, W., Protheroe, R.J. 1997, PRL, 79, 2616 
Bednarek, W., Protheroe, R.J. 2001, Astropart.Phys., in press 
Berezinsky, V.S., Prilutsky, O.F. 1978, A&A, 66, 325 
Bertou, X., Billoir, P., Deligny, O., Lachaud, C, Letessier- 

Selvon, A. 2001, Astropart.Phys., in press 
Blasi, P., Epstein, R.I., Olinto, A.V. 2000, ApJ, 533, 123 (BEO) 
Burrows, A., Lattimer, J.M. 1986, ApJ, 307, 178 
De Goubeia Dal Pino, E.M., Lazarian, A. 2000, ApJ, 536, L31 
Fichtel, C.E., Trombka, J.I. 1981, Gamma Ray Astrophysics, 

(NASA SP-453, Washington, DC), 192 
Giller, M., Lipski M. 2001, Proc. 28th ICRC (Hamburg), 

p. 2092 

Kobayashi, C, Tsujimoto, T., Nomoto, K. 2000, ApJ, 539, 26 
Learned, J.G., Mannheim, K. 2000, Ann.Rev.Nucl.Part.Sci., 
50, 679 

Lindblom, L., Owen, B.J., Morsink, S.M. 1998, PRL, 80, 4843 
Mannheim, K., Protheroe, R.J., Rachen, J. P. 2001, Phys.Rev., 

D63, 23003 
Narayan, R. 1987, Ap J, 319, 162 



Protheroe, R.J. 1998, astro-ph/9809144 



Protheroe, R.J., Bednarek, W., Luo, Q. 1998, Astropart. Phys., 
9, 1 

Spiering, C. 2000, Nucl.Phys.Suppl., 91, 445 
Waxman, E., Bahcall, J. 1999, Phys.Rev. D, 59, 023002 
Wheeler, J.C., Yi, I., Hoflich, P., Wang L. 2000, ApJ, 537, 810 
Woolsey, S.E., Langer, N., Weaver, T.A. 1993, ApJ, 411, 823 



Acknowledgements. I thank the referee for comments and sug- 
gestions. This work is supported by the Polish KBN grant No. 
5P03D02521. 



